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The violations of the fluctuation- dissipation theorem are analyzed for a trap model with 
a gausssian density of states. In this model, the system reaches thermal equilibrium for 
long times after a quench to any finite temperature and therefore all aging effect are of 
a transient nature. For not too long times after the quench it is found that the so-called 
fluctuation- dissipation ratio tends to a non-trivial limit, thus inicating the possibility for the 
definition of a time scale dependent effective temperature. However, different definitions of 
the effective temperature yield distinct results. In particular plots of the integrated response 
versus the correlation function strongly depend on the way they are constructed. Also 
the definition of effective temperatures in the frequency domain is not unique for the model 
considered. This may have some implications for the interpretation of results from computer 
simulations and experimental determinations of effective temperatures. 



PACS numbers: 05.40.-a, 64.70.Pf, 61.20.Lc 



I. Introduction 

The out-of-equilibrium dynamics of glasses and disordered systems has been studied quite 
intensively in the last decade, see e.g.pp. One reason for the interest in this area lies in the 
fact that glassy systems usually do not reach thermal equilibrium within experimentally 
accessible times when cooled to low temperatures. In such a situation the question as 
to whether there is a possibility for characterizing the non-equilibrium state of the sys- 
tem naturally arises. In particular, it would be extremely helpful if some of the tools of 
statistical physics could be extended to non-equilibrium situations. One of the efforts of 
quantifying the deviations from thermal equilibrium is concerned with extensions of the 
fluctuation-dissipation theorem (FDT) relating the response of the system to the two-time 
correlation function. In equilibrium, this relation is determined by the thermodynamic 
temperature. For non-equilibrium systems, the deviations from FDT can in some cases be 
used for the introduction of a time scale dependent effective temperature Teg [2]. A vari- 
ety of model-calculations and numerical simulations have been performed to analyze the 
behavior of Tes obtained this way, for reviews see [Sill]. Recently, also some experimental 
investigations of the violations of the FDT have been performed, see e.g.[5j. One typically 
finds a behavior in which Tee exceeds the bath temperature by a certain amount. This 
makes sense as the typical preparation of the system is given by a quench from a high 
temperature into a glassy phase. 
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There are, however, still some open questions regarding the usefulness of the inter- 
pretation of the effective temperatures obtained from FDT-violations. For instance, for a 
meaningful definition of a temperature, Tes should be independent of the dynamical vari- 
able used for its calculation. In case that so-called neutral variables [U] are used. Teg often 
is found to be independent of the variable in some long-time limit. Another question is 
concerned with the relaxation of Tgfj to the bath temperature when the system reaches 
thermal equilibrium on some time scale. This point has been investigated in some recent 
papers [3 IHl E] with the result that some care has to be taken in the determination of Tes- 

In the present paper, I will consider another model in which the system reaches equi- 
librium for long times after a quench, namely a trap model with a Gaussian distribution 
of trap energies|Tn]. The same model has been introduced earlier as an approximation to 
a random walk model^I] and was termed 'energy master equation' then. It was shown in 
ref.pi] that this model is able to reproduce some features that are observed in the relax- 
ation of supercooled liquids and glasses. Furthermore, recently it was demonstrated that 
some results of simulations on models of supercooled liquids can (at least partly) be inter- 
preted in terms of a Gaussian trap model[T21 UHl E]. Some features of the aging behavior 
of the populations of the traps have already been discussed in ref. [7] . There it was shown 
that after a quench from an (infinitely) high temperature the distribution of populations 
as a function of the time elapsed after the quench first narrows and then broadens again, 
in qualitative accord to some observations in computer simulations JS]. 

In the following Section, I will briefly recall the derivation of the fluctuation-dissipation 
relations for trap models^J El CHI- Section III the results of model calculations are 
presented and discussed. The paper closes with some conclusions in Section IV. 

II. Fluctuation-dissipation relations for trap models 

The derivation of the fluctuation-dissipation relations for arbitrary continuous-time Markov 
processes has been given in ref.pH]- Usually one assumes a situation in which the system 
under consideration is quenched from a high temperature Tq to a low temperature T in 
the beginning of the experimental protocol. After a time t^, called the waiting time, has 
elapsed after quench, some dynamical quantities are monitored. In the present context, 
the correlation function of some dynamical variable M{t), 



and the response to a field applied at time t^, H{t) = H6(t — tw), conjugate to M, 



are of particular importance. In the above expressions, the values of e are meant to 
represent the trap energies. M{e) is the value the dynamical variable M acquires in 






(2) 
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the trap with energy e and G{e,t\e' ,t^) is the conditional probabihty to find the sys- 
tem in trap 'e' at time t, provided it was in trap 'e" at time ty^. The populations 
p{e',tu]) evolve from the initial populations p{e,t = 0), typically equilibrium populations 
at the starting temperature Tq, via p(e, t^) = Jde'G{e,tw\e' ,0)p{e' ,0). For a stationary 
Markov process, the conditional probability G{e,t\e' ,tw) is time-translational invariant, 
G{e, t\e', tin) = G{e, t — t^|e', 0) = G{e, t — tw\e'), and obeys a master equation P^: 

^G(e, t\t') = - Jde"W{e"\e)G{e, t\e') + Jde"W{e\e")G{e" , t\e') (3) 

with W^(e|e') denoting the rates for a transition from e' to e. In order to calculate the 
linear response to a field one has to fix the dependence of the transition rates on the field 
H. Following RitortpTC], I use the following form of multiplicatively perturbed transition 
rates: 

iyW(e|e') = iy(e|e')e^^[^^W-/^^(^')l (4) 

where 7 and fi are arbitrary parameters and (3 = 1/T (with the Boltzmann constant set 
to unity). If /i + 7 = 1 holds additionally, then the rates W^^\e.\e') obey detailed balance 
also in the presence of the field, provided the Vr(e|e') do so. In all following calculations 
the unperturbed transition rates will be chosen according to: 

W{e\e') = r]{e)K{e') with K{e') = K^e^^'' (5) 

as is usual for the trap modelfOl El 120] ■ Note, however, that also other choices have 
been considered, see e.g. As already mentioned above, I will only consider the trap 



model with a Gaussian density of states, 77(e) = J- e in this case, the system 



/2-KCJ 

reaches equilibrium for long times regardless of the initial conditions. The equilibrium 
populations are given by p'^^(e) = limj^oo G(e, t|eo, to) = ^/^^ e"^^"^-*^/^^""^^ with the tem- 
perature dependent mean energy e = — /Scx^ITT]. Note that the variance a is temperature 
independent. 

The details of the calculation of the linear response according to eq.Q have been 
presented in refs.jTHl El UHl- As a result, it is found that for mean-field trap models with 
a dynamical variable M that obeys a distribution with zero mean, (M) = 0, and unit 
variance, (M^) = l, the fiuctuation-dissipation relation can be written as: 



Here, the correlation function 



9n(t,t^) 9n(t,t, 



7 — Tz: 



dt,,, dt 



(6) 



U{t,t^) = dee-<'^^'-'-^p{e,t^) (7) 



gives the probability that the system has not left the trap occupied at in the following 
time interval (t — t.„)[TU]. 
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In equilibrium, all quantities are time-translational invariant and one has 

rfneq(t) 



i?eqW = -/5(7 + /i)- 



dt 



which for /i = 1 — 7 is just the well known FDT. Note that for the present choice of the 
transition rates, eq.©, the system always reaches equilibrium for ty^ —>■ 00. Therefore, all 
aging effects are of a transient nature. This behavior is similar to what one expects for 
structural glasses cooled not too low below the calorimetric glass transition temperature. 
In the present paper, the quenches are performed in the following way: In the beginning, 
the system is prepared in thermal equilibrium at the initial temperature Tq, i.e. P2?(e) = 
1 g-(e-»;o) /(2cr ) -^ji;]^ g _ _g fjj^g following time evolution is calculated at the 



working temperature T, pr(e, ^«)) = /c^e'Crle, tiu|e', O)p^(e) where Grle, 0) is the 
solution of cq.® with the transition rates (0) evaluated for /? = 1/T. 

The violations of the FDT in out-of-equilibrium situations can be characterized by the 
so-called fluctuation-dissipation ratio (FDR): 

which also allows the definition of a time scale dependent effective temperature Teg (t, t^) = 
T/X{t,tw), the detailed behavior of which has been discussed for a variety of modelsjS]. 

Another way to define an effective temperature has been proposed and applied in ex- 
perimental determinations of FDT- violations [21 E]: 

T..M = '^^^j^ (10) 

relating the Fouriertransform of the correlation function (spectral density) to the dissi- 
pative part of the susceptibility. Of course, one would assume the various definitions to 
yield identical results. In the next Section, different definitions will be compared and the 
differences will be discussed for the Gaussian trap model. 



III. Results and Discussion 

I will start with a brief discussion of the equilibrium properties of the correlation function 
function n(t, t^). As already mentioned above, the system reaches equilibrium for long 
times after a quench has been performed, i.e. in the limit 00. According to eq.(|7j). 

in this limit one has for t = {t — ty^) < 00: 

neq(r) = |dee-''(^V^(e) (11) 

This function is plotted in Fig.l as a function of r/roq for various temperatures, where r^q 
is the relaxation time determined as the 1/e decay time of neq(r). It is evident that the 
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decay becomes broader with decreasing temperature. This means that in this model time- 
temperature superposition does not hold in equilibrium. Without showing the results here, 
I only mention that the values for the relaxation time Teq and the stretching parameter 
Peq in a Kohlrausch fit to Ileq, neq(r) ~ exp — (T/req)^°i , roughly follow the expressions 
given by Monthus and BouchaudpH]. In ( Teq) ~ (1/7")^ and Peq ~ [1 + (o./T)]^^^'^ where a 
is a constant. 

In order to see in which time range aging effects can be observed in n(t,t^) for finite 
t^, in Fig. 2a I show n(t^ + r, ty^) as a function of the measuring time r for a quench from 
infinite temperature, /?o = 0, to various working temperatures. It is seen that the decay 
times T{tw) span a much broader range at low temperatures than at high temperatures. 
A closer look additionally reveals the fact that the decay is broader for small and large 
waiting times than for intermediate tw The origin of this behavior lies in the waiting 
time dependent width of the distribution of the population, p(e, t^), cf. the more detailed 
discussion on this point in ref.jZj. 

The time window in which significant aging effects are to be expected can most easily 
be quantified by considering the mean relaxation time. 



which does not depend on the width of the actual decay. The mean relaxation times for a 
quench from Tq = oo to various temperatures are shown in Fig.2b. The difference of the 
mean relaxation time directly after the quench, (tq) = (t(0)) = kJ^^ exp (/3^(T^/2) and the 
equilibrium mean relaxation time (rgq) = (r(oo)) = exp (3/5^(7^/2), shown in the inset 
in Fig. 2b, gives a measure for the 'aging time window' at a given temperature. 

As is exemplified in Fig. 3a, the overall time window for visible aging effects diminishes 
not only with increasing final working temperature, but also with lower initial temperature 
Tq. It is evident that the effect of lowering Tq is similar to measuring the correlation 
at a higher working temperature for a given Tq. In Fig. 3b the relaxation time 'r(t^) 
and the stretching parameter I3{tw) as obtained from Kohlrausch fits to n(t^ + T,tw), 
n(t^ + T, t^) ~ exp — (T/r(tu,))^(*™) , are shown. One can see that the overall spread in 
relaxation time diminishes with decreasing Tq and that r(t^) approaches Teq at the working 
temperature T on the time scale of Teq itself. This is because there is no other time scale 
present in the model. The stretching parameter (3{tw) first increases and then decreases 
again as a function of the waiting time, similar to what has been found earlier in a free- 
energy model for the primary relaxation of viscous liquids jlH]. Note that in the limit of 
short and long the /3(ti„)-values, /?(0) = /3(cxd), are all the same because in the present 
model the width a of rj{e) and p'^'^{e) are independent of temperature. 

Next, we turn to a discussion of the FDR, eq.®, which is calculated from the 
fiuctuation-dissipation relation, eq.®. For long waiting times one finds 

t^) = 7 + /i for ty, > Teq (12) 

as expected because the system is in equilibrium. Furthermore, also the equal-time value 
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yields the same result, 

X{t,t)= lim X(t,t^)=7 + /i (13) 

This can easily be seen from the fact that for small time-differences one has dt^Il{t,tw) = 
—dtll{t,tw). For finite and long times, however, the FDR reaches a non-trivial limiting 
value 

Xoo(t^) = ^ lim X(t,t^)=7 (14) 

independent of Tq > T. This result is very similar to the one obtained in ref. 0|. It is 
interesting to compare the result with the corresponding expression for the trap model 
with an exponential density of states, where it is found that the long time limit of Xoo{tw) 
tends towards 7, X^o = lim^ _^ X^ft,.)P^. Thus, apparently, the model with a gaussian 
density of states resembles some features of the model with an exponential one for interme- 
diate waiting times. If interpreted in terms of a time scale dependent effective temperature, 
T/Xoo{tw) = 1/7, one finds that it can be higher (7 < 1) or smaller (7 > 1) than the bath 
temperature T. It should be mentioned, however, that one typically will have 7 < 1. For 
instance, in a 'force model', one would assume the bias induced by the external field to be 
symmetric and therefore 7 = /i = 1/2. 

The detailed behavior of the FDR is shown in Fig.4, where I plotted X{tw + T, t^) versus 
the measuring time r for various and 7 = for a quench from Tq = 00 to T = O.Scr. It 
is seen that X{tw + T^t^) starts from the short-time value X(t^,t^^,) = 1 and decays to its 
long-time limit, eq.lfT^. This decay takes place on the time scale (r/t^) ~ 1 as is to be 
expected for a model with a single time scale. For very long t^, the decay of X will not be 
observed due to the long measuring times required. If one starts with lower Tq the main 
difference is that the decay from X = 1 to X = 7 for a given takes place at somewhat 
later measuring times r. This is because for smaller Tq the system is closer to equilibrium 
for a given t^. Thus, one roughly can compare X for the lower Tq and a given t^, with the 
corresponding one for the higher Tq but a longer t^. 

Instead of computing X{t,tyj) directly, the FDR or the effective temperature is usually 
obtained from a so-called fluctuation-dissipation (FD) plot, i.e. a plot of the integrated 
response x{t-,tw) = JljlsR{t,s) versus the correlation function n(t, t^„) [21 E]- The FDR 
then is extracted as the slope in such a FD-plot. It has been pointed out earlier that such 
plots have to be constructed with some care in the general casejH]. In some situations 
a plot of x(tw + Tytw) vs. Il(tw + T,tyj) with r as the curve parameter and fixed tu, 
does not necessarily give the correct FDR. This is only ensured if is used as the curve 
parameter due to the fact that according to the definition of the integrated response one 
has R{t,t^) = -^2c^_ The reason is that in general the different derivatives dt^Tl{t^tw) 
and djl{t,t^) cannot be interchanged. For some models with transient aging behavior, it 
has indeed been found that the former construction yields wrong results jHl [HI. Also for the 
gaussian trap model, FD-plots with r as the curve parameter yield wrong results for the 
FDR. This fact is exemplified in Fig. 5a, where such FD-plots for different temperatures and 
tu, = lO^^^Tcq are shown. The initial temperature was chosen as Tq = 00. The inset shows 
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the temperature dependence of the slopes extracted from these plots via a linear regression 
in the interval 0.3 < 11 < 0.7. Note that all slopes are larger than Xoo(t«,) = 7 = 0.5. The 
situation is very similar to the one in a free-energy model for glassy relaxation discussed 
in detail in ref.jU]. 

In Fig. 5b, a FD-plot with as the curve parameter and various values of the time 
t is shown for 7 = (upper panel) and 7 = 1/2 (lower panel). It is evident that in 
this case the limiting slopes coincide with Xoo(t«,), eq. lfT^ . However, it is also clear that 
the curves change continuously until this limiting slope is reached. In a strict sense this 
would mean that a thermometer would measure an increasing effective temperature as a 
function of time even though in the model there is only a single time scale. Therefore, 
similar to the case of the trap model with an exponential density of states, it appears 
that only the limiting slope Xoo{tw) may be a useful candidate for the definition of an 
effective temperature jn|. Additionally, the detailed behavior of the curves also depends on 
the initial temperature Tq, as can be seen from Fig.5c which shows FD-plots for various Tq 
and a fixed time t = 10~^rcq for 7 = 0, ;U = 1. The limiting value for the slope is the same, 
as already pointed out above. 

From the above discussion it appears that the quantity T/Xoo(tw) may serve as a 
t^-dependent effective temperature in the gaussian trap model. As already noted in the 
previous section, there are other possible definitions of effective temperatures, in partic- 
ular the one given in eq. lfTIH) . Tes{uj,tw) = ^^^t^yj- Here, the natural way to define the 
Fouriertransform is with respect to the measuring time 



C{u}, t^) = Re / dTU{t^ + r, t^)e*^^ 
Jo 

X"iuj, t^) = Im / drRit^ + r, t^)e'^^ (15) 

JO 

which is similar to the definition used in most experiments. It can be shown analytically 
that Tcff(u;,t^) defined this way has the following limiting behavior: 

T 

lim Tcs{uj,tnj) = lim \im Tcs{uj,tuj) = — ■ — (16) 

t^^oo w^oot^^O /i + 7 

as expected. However, the low-frequency hmit is given by: 

T 

lira }im T^s{uJ,t^) = — ■ ^^rj^ (17) 

which is different from T / X^{tyj). For /i = 1 and 7 = no deviations from the equilibrium 
FDT can be observed at all. When considered as a function of frequency, Tefj(a;,t^) starts 
from the low-frequency limit T / (/i + 76^^°"/^) ) and roughly at ut^, ~ 1 smoothly crosses 
over to the high-frequency limit T/(/i + 7). The main difference to T/X is that the low- 
frequency limit usually is much smaller than the bath temperature T. An obvious reason for 
this discrepancy lies in the definition of the Fouriertransform. The definition of the FDR 
involves only the partial derivative 9j^H(t, t^^,), whereas ujC{(jj,tyj) according to eq. (fT3|) 
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corresponds to the the Fouriertransform of dtll(t,tw)- Thus, the situation appears to be 
very similar to the one met when considering the different ways of constructing FD-plots 
discussed above. There is, however, one difference. For the present model it is not at all 
clear that a definition of an effective temperature via the Fouriertransform of n(t, t^) and 
R{t, tyj) is related to the FDR in a direct manner. This is due to the fact that one has to 
consider the Fouriertransform of R{t,t^) = pX{t,tw)dt^Il(t,tw) which in general will give 
a convolution instead of a simple product. 

It is tempting to use an alternative definition of the Fouriertransform [2], C{uj,t) = 
Re /Q(isn(t, s)e~*'^(*~'') and x"(uj,t) = Im jQdsR{t, s)e~'^'^^^~^\ for which one can see that 
uC{u, t) is related to the Fouriertransform of dt^Il{t, t^). One then has Teff (cu, t) = '^^(^f) • 
However, this definition implies that one considers the situation cut 3> 1, in order for the 
oscillatory part of the functions involved to be negligible. As one can show analytically 
that Tefr(cc;,t oo) = Tefr('^ ^ oo,t) = T/{j + /i) this condition means that the relevant 
small frequency-limit is hard to reach in practice, in particular because the relaxation takes 
place on the scale cut ~ 1. 



IV. Conclusions 

In the present paper I considered the aging behavior of the correlation function and the 
linear response in a gaussian trap model. As in this model equilibrium is reached for 
long waiting times after a quench (t^ oo), all aging effects are of a transient nature. 
The time window in which aging effects can be observed expands with decreasing working 
temperature T if a quench from a fixed high temperature Tq is considered. If T is kept 
fixed and Tq is increased the 'aging time window' also increases. 

If the waiting time is chosen to be longer than the equilibrium relaxation time (1/e 
decay time), ^ Teq, the two-time quantities like the correlation function n(t, t^) and 
the response function R{t, tw) approach their equilibrium values and only depend on the 
time-difference {t — tuj). In this situation the FDT is obeyed. For short waiting times, on 
the other hand, one has an explicit dependence of n(t,t^) and R{t,tw) on both times and 
the functions no longer are time-translational invariant. In this situation strong violations 
of the FDT are observed. In particular, for < Teq the FDR X{t,tw) tends to a non- 
trivial long-time limit Xoo(t«,) = 7, cf. eq.(fT^. independent of Tq > T. The constant 7 is 
determined by the coupling of the system's dynamics to an external field via the transition 
rates W^^\e'\e), cf. eq.(IH). The only possible temperature dependence of Xoo{tw) stems 
from a T-dependence of 7. The value for Xoc{tw) coincides with that found for the trap 
model with an exponential density of states JHl; albeit for long waiting times. 

The FDR can be used for the definition of a time scale dependent effective temperature, 
T/X{t,tw). Usually, the effective temperature is extracted from FD-plots, which can be 
constructed in different ways. As has been pointed outjO], the correct way to construct FD- 
plots is to use tw as the curve parameter for fixed t. If FD-plots are constructed this way, 
the limiting slope coincides with Xoo(t«,). The fact that these plots are curved indicates 
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that not the function X{t, t^,) as a whole but only the limiting value Xoo(t«,) may be used 
for a meaningful definition of an effective temperature. If alternatively FD-plots with the 
measuring time as curve parameter for different t^, < Teq are considered one finds that 
one can define a slope in these plots which, however, does not coincide with and 
furthermore show some artificial temperature dependence. A similar problem occurs when 
the definition of an effective temperature in the frequency domain is considered. The usual 
definition of the Fouriertransform, eq. ()15|l . yields an effective temperature that is distinct 
from what one finds from the FDR. In this case, similar as for the FD-plots, one reason is 
that one cannot interchange the derivatives with respect to the earlier and the later time 
in the correlation function. 

In the recent past FD-plots have been constructed from computer simulations of model 
supercooled liquids|U For the model of a fragile liquid one finds that the FDT is 

obeyed for large values of the correlation and then the slope changes to a smaller value 
indicating an effective temperature larger than the bath temperature jH For a strong 
liquid, however, one finds an effective temperature smaller than the bath temperature |25j . 
These plots have been constructed with the measuring time as the curve parameter and 
fixed waiting times. 

Given the fact that the gaussian trap model apparently is able to capture some impor- 
tant features of the dynamics observed in computer simulations on supercooled liquids, it 
is tempting to speculate that FD-plots with the waiting time as the curve parameter might 
look different to the ones published so far. Therefore, it would be interesting to calculate 
FD-plots from simulation data the correct way and compare the results to the previous 
findings. 
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Figure captions 



Fig.l : neq(T) versus r for various temperatures. 

Fig.2 : a: Il{tw + r.tw) versus measuring time r, scaled to the equilibrium relaxation 
time Teq (1/e decay time) for a quench from Tq = cxo. 

b: {T{tw)) versus tw/{Teq) showing the time window in which aging effects are to 
be expected. The inset shows the mean relaxation times (tq) = (t(0)) and (req) = 
(t(oo)) as a function of inverse temperature. 

Fig.3 : a: Il{ty^ + T,tyj) versus r/rcq for quenches from various Tq = 0, a, 0.5a and 

^ogioit^/r^q) = -15, -5, -3, -1, 1, 5. 

b: T{tw) (upper panel) and P{tw) (lower panel) versus t^/Teq as obtained from fits to 



a Kohlrausch function, exp —(T/T{tyj))^^*' 



Fig.4 : The FDR X(tw + r, t^) versus r/req for a quench from Tq = oo to T = 0.3a and 
7 = 0, /X = 1 for various waiting times t^. 

Fig.5 : a: Plot of Tx{tw + 't, tw) versus n(t„,, + r, t^) (FD-plot) with the measuring time 
r as the curve parameter and fixed waiting time, = lO^^^Tcq, for a quench from 
To = oo to T = 0.2, 0.3, 0.4, 0.5, 0.6, 1.0, 3.0a from bottom to top. Here, 7 = /x = 1/2. 
The inset shows the slopes extracted from these plots via a linear regression for the 
data in the interval 0.3 < 11 < 0.7. 

b: Plot of Tx{t, tyj) versus U{t, tw) (FD-plot) for a quench from Tq = 00 to T = 0.25a 
and times \ogiQ{t / Teq) — —6, —4, —2,0, 6 from bottom to top. In the upper panel 
7 = 0, /i = 1 has been chosen, i.e. Xco{tw) — and in the lower panel j — /i — 1/2 

(Xoo(t«,) = 1/2, dotted line). 

c: Tx(t,t^) versus n(t, t^u) for a quench from Tq = 00, l.Ocr, 0.5a, and 0.4cr to 
T = 0.25a and t/r^q = 10-^ 7 = 0, = 1. 
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G. Diezemann, Fig.1 
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G. Diezemann, Fig. 2a 
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G. Diezemann, Fig. 2b 
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G. Diezemann, Fig. 3a 
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G. Diezemann, Fig. 3b 
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G.Diezemann, Fig. 4 
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G. Diezemann, Fig. 5a 
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G. Diezemann, Fig. 5b 
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G. Diezemann, Fig. 5c 
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